Do quantum effects hold together DNA condensates? 
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Abstract 

The classical electrostatic interaction between DNA molecules in 
water in the presence of counterions is reconsidered and we propose 
it is governed by a modified Poisson-Boltzmann equation. Quantum 
fluctuations are then studied and shown to lead to a vacuum interaction 
that is numerically computed for several configurations of many DNA 
strands and found to be strongly many-body. This Casimir vacuum 
interaction can be the "glue" holding together DNA molecules into 
aggregates. 
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The theoretical study of the formation and stability of DNA aggregates 
achieved an important understanding of these phenomena that, nonetheless, 
still appear mysterious in many respects pp. A DNA molecule in aqueous 
solution ionizes and gives rise to a highly charged anion that, in the presence 
of cations, binds them by the Oosawa-Manning (OM) condensation [2, 3j (for 
a review see, e.g., [TJ HI [5]). When about 90 per cent of the DNA negative 
charge is screened and when the cations have a specific valency +k (usually 
k = 3 or k = 4) the DNA strands collapse to form rod-like, spheroidal and 
toroidal aggregates [6]. 

There are two approaches followed to understand these features: In one 
the detailed charge distribution of the DNA is used to calculate the elec- 
trostatic interaction between two DNA strands using a linearized Poisson- 



Boltzmann (PB) (or Debye-Hiickel (DH)) equation [7] (see also [8]). Within 
this approach the attraction between two like-sign charged but suitably ori- 
ented DNA helical strands and the specificity of cations driving the attrac- 
tion can be predicted. However, it was also realized that an assembly of 
strands in an hexagonal array is a frustrated-charge system [9]. Thus the 
reasons why within this approach the DNA aggregates should form and be 
stable is not clear, unless other forces are invoked. 

In alternative approaches the surface of a single DNA strand is treated as 
a two-dimensional complex system and statistical mechanical arguments lead 
to the counterion-mediated attraction between DNA strands [101 [TTj \12\ [T3] . 
The key idea there is that the condensation is triggered by local correlations 
and thermal fluctuations not present in the mean-field PB approach. The 
electrostatic interaction again plays an essential role and an assembly of 
DNA strands forming an hexagonal bundle is found to be a frustrated-charge 
system |14| . Here, again, the formation and stability of aggregates is not 
clear if other forces are not present. 

In this letter we propose that quantum vacuum fluctuations are likely to 
be the solution to some of these puzzles. In particular we propose that the 
Casimir interaction due to these fluctuations is the "glue" that holds DNA 
aggregates together. 

The PB equation we need to consider is 

n^.^ Air m . 87T , fke<f>(x)\ , . 

V 2 $(f) = -— p(g,T) = — kenosmhi^-^J , (1) 

where <3?(x) is the electrostatic potential due to the DNA strand (seen as a 
negatively charged rod immersed in water at room temperature, T ~ 300K, 
with dissolved salt whose ions have valency z = ±k, with k = 1,2,...) and to 
the ions, the medium has dielectric constant e, and the charge distribution 
of the composite system DNA-salt is p(x,T) = pbna(x, T) + ke(n + (x,T) — 
T)), with n±(x,T) the concentration (density) of ions following a 
Boltzmann distribution. Usually pdna is not included into the PB equa- 
tion ([1]). We demand, instead, that it obeys a Boltzmann distribution law 
as for the ions [T5] 

PdnaO?, T) = -T»DNA(x)|g| exp ( 2 ) 

where q < is the charge of the DNA strand with n^ NA (x) = J2iLi l/ i(zi)6 2 (x±— 
li) and N the number if strands. This charge density function also defines 
our approximations: we model the DNA strands as infinite lines all parallel 
to the z-axis and located at li in the x — y plane with the coefficients V{{zi) 
carrying information on the charge structure of the DNA strand. We further 
simplify our model by taking Vi(zi) = v = constant, Vi = 1, N. 

Our concern is to study the interaction among DNA strands after at 
least 90 per cent of the negative charge has been screened via the OM 
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Figure 1: The configuration used for 19 DNA strands where x is the distance 
between nearest neighbors. 

condensation [2], [3], as this is the reported critical value for collapse. It is 
then reasonable to consider <E> small. Thus, expanding the exponentials till 
the first order we obtain 



N 

' ' Hx) = J , (3) 



i=l 



which is a modified DH equation. Here = k 2 K 2 , with k 1 = (efcsT/(87re 2 no)) 1 / 2 
the Debye screening length, A = 47r^|g| 2 /e/csT, and J = — (l/e)4>jr\q\u Yli=i ^ 2 \x±- 
li). 

We now consider small time-dependent fluctuations, 3>(x) — > &{x) + 
4>(x,t), where <3? satisfies Eq. (|3|) that descends from the action 



\ 1=1 



x ± -li)]$ + J$) , (4) 



where we use units h = c = 1, with c the velocity of light in the medium and, 
for the sake of clarity, we included an integration over time Jq dt even though 
the functions are time-independent. We then demand that to the fluctuation 
field (j> as well is associated an action that is a suitable modification of 
namely 

A{4>) = |^^^-a 2 -9 2 -Vi + ^ 2 + A^( 2 )(x ± -0^. (5) 

Note that in A(4>) the term with the coupling to the "external current" J is 
zero because / d A xJ(j) = J d?xJ f£ dtrf* = as required for fluctuating fields. 
The way to consider the effects of <f> is to average these fluctuations out to 
obtain an effective action A e s(§). This is done by considering the generating 
functional Z[<&,4>] = J[D<$>]e iA W f[D<t>] e iA ^ = f[D$] e -(-4(*)+corrections) 
where we Wick rotate on the time direction t — > it, and identify A e g(&) = 
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Figure 2: Energy of interaction of two DNA strands. The lower (upper) 
curve corresponds to a = 2 (a = 1). Distances are measured in units of 
~ O(10) A, while energy units are estimated to be ©(lO" 1 ) - 0(1) eV. 



A($)+ corrections. The result is [E] Aff(^) = A($) + Sr, with £ = 
3 /- ~ fl Jo + °° dEp{E)^E + p 2 , where is the density of states. 

This energy is of the form £ = (1/2) w > i- e - it is the zero point Casimir 
energy of the system and it can be determined [TS], [TB] (see also [T7]) and 
is 

fir r°° 
£ = ^ d Eln 

where we reintroduced H and c, Kq{x) is the modified Bessel function of 
the second kind of order zero, \x is the mass (inverse length) scale parameter 
introduced earlier, M is a further mass (inverse length) scale parameter that 
satisfies M < jx and lij = \h — lj\ are the relative distances between DNA 
strands. Note that A, the other parameter of the theory, is first renormalized 
A — > Xr and then adsorbed into the definition of M — > Mexp(27r/A^) = M. 

Our strategy is to study the energy of configurations of DNA strands that 
capture as much as possible the symmetry of arrangements encountered in 
real cases [B]. Hence, after having learned on the two-strand interaction, we 
focus on many-body interactions where the strands are sitting at the sites 
of hexagonal lattices^ (see Fig. [T]for the case of 19 strands) and perform a 
careful analysis. 

To render the expression © suitable for such a study we first numerically 
perform the integral over E and then plot the resulting expression £{x) 
where the relative distances 1™ = CijX are expressed in terms of the basic 
lattice distance x, the numerical coefficients Cij take the symmetry of the 

1 We also present here results for four DNA strands sitting at the vertices of a rhombus. 
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Figure 3: Interaction energy for the hexagonal lattice 19-strand configura- 
tion of Fig. [TJwith a = 1 (upper curve) and a = 2 (lower curve). The units 
are ~ 0(lO)Afor the distances and O(10 _1 ) — 0(1) eV for the energy. 



given arrangement into account, the distances are measured in units of /i ~ 
O(10) A, the other scale M is constrained to be positive and less than //(= 1) 
and we write it as < M = e~ l l a < 1. In this fashion the range of £ scales 
with a and we present here results for a = 1 and a = 2. To estimate the order 
of magnitude of the energy in the range of non- vanishing interaction we write 
the integral in Eq.© in a dimensionless fashion introducing the scale \x. This 
way we obtain that the energy scale factor in front is hcfi ~ 6 x 10 2 eV, where 
c is taken to be the speed of light in water, c ~ 2xl0 8 ms _1 . Considering that 
the numerical values of the integrals we obtain are 0(1CT 3 ), the estimate of 
the magnitude of the effect in the range is £ ~ 0{1Q~ 1 ) — 0(1) eV, that is 
between one and two orders of magnitude stronger than the thermal energy 
that at room temperature is fc^T ~ 2 x 1CP 2 eV. These estimates clearly 
indicate that the quantum relativistic effect we are considering is important 
for explaining the attraction of DNA strands. 

The two-strand interaction energy is shown in Fig. It is clearly at- 
tractive and finite-range. Similar attractive behaviors for the two-body in- 
teraction have been found in various models [TT] . What we observe here 
is that in those models it is not clear why the interaction still needs to be 
attractive for more than two strands and why the aggregates are stable. For 
the Casimir energy we are considering here this is indeed the case, since this 
attraction mechanism does not suffer of any frustration. 

To establish whether the magnitude and range of this attractive energy is 
indeed relevant for the case of DNA aggregates we need to move to the many- 
body case. We have computed the energy for several interacting strands 
having various configurations. We present in Fig|3]the results for 19 strands 
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Figure 4: The lower curve is the interaction energy for the four-strand- 
rhombic configuration (many-body). The upper curve is what is obtained 
by summing-up the 6 two-body interactions. In both cases a = 2. The units 
arc ~ O(10) A for the dist ances and O(10 x ) — 0(1) eV for the energy. 



arranged as in Fig.[TJ More configurations are discussed in [IS]. Comparing 
these plots with that of the two strands interaction we clearly see that the 
attraction becomes stronger and acts on a larger range when the number of 
strands increases. The range of attraction can be adjusted by fixing a to fit 
the typical distances reached within the aggregates, that for the hexagonally 
packed toroidal condensates ranges between [6] 18 A and 28 A, values clearly 
compatible with the range we obtain here. 

In real cases it is always several DNA strands that interact, the two 
strands being only an idealization. Thus the fact that for 19 strands we find 
that (for a = 2) the range of the force is in agreement with the typical values 
reported for DNA aggregates [6] we take it as an indication of the validity of 
our hypothesis that the quantum Casimir energy holds together the aggre- 
gates. Furthermore, this force is many-body in nature and the many-body 
effects are big, another reason for taking the two-body interaction only as an 
indication of the real phenomenon. That the many body effects are strong 
we proved in our numerical calculations where we compared the iV-body en- 
ergy of Eq. ([6]) with that obtained by summing up (N/2)(N — 1) two-body 
interactions. The results for four and seven strands are shown in Fig. 0] and 
Fig. [5j respectively, and they indicate that the effect grows with N. 

The singularity in the energy at a value x such that the determinant 
function becomes zero, is the indication of the limit of validity of our ap- 
proximations, x can be evaluated for the various cases by plotting the de- 
terminant [15] . That singularity means that if only the Casimir force were 
present the strands would collapse to zero separation, an instance that does 
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Figure 5: The lower curve is the interaction energy for the 7-strand con- 
figuration, i.e. for 7 strands sitting at the vertices and at the certer of 
a regular hexagon (many-body). The upper curve is what is obtained by 
summing-up the 21 two-body interactions. In both cases a = 2. The units 
are l//z ~ 0(1O)A for the distances and 0(1O _1 ) — 0(1) eV for the energy. 

not occur in the real case because of nonlinear corrections to the Casimir 
force itself at every short distance and because of other forces (not consid- 
ered here) such as the electrostatic force that for more than two strands will 
give a net repulsive effect. Another important factor at such short distances 
is of course the finite size of DNA strands that have a transverse length 
(radius of the cylinder) of 10A. 

The main result we present here is probably the demonstration that 
quantum relativistic effects can be responsible for the collapse of DNA strands 
into aggregates (after the OM condensation has taken place) and for holding 
them together into stable condensates. The time-dependent fluctuations of 
the electric field that we studied are quantum in nature, propagate at the 
speed of light in the medium thus give rise to a Casimir force that is attrac- 
tive and short range for the two-body case and is many-body, the departures 
from the "sum of two-bodies" being important and growing with the num- 
ber of strands. The magnitude and range of the interaction is such that it 
could explain the formation and stability of DNA aggregates, as a prelimi- 
nary comparison with reported data shows. We focused our attention on the 
difficult problem of computing such interaction for several DNA strands and 
were able to overcome the analytical challenges with numerical calculations 
performed for a variety of cases, most of which with hexagonal symmetry of 
arrangement, the typical situation reported in experiments. Although our 
model is a simple one and to have the full picture it needs to be completed 
(by introducing other forces such as the hydrophobic, hydrophilic and elec- 
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trostatic forces, and finite size effects, as well as by considering nonlinear 
corrections to the Casimir effect itself at very short distances x) in the light 
of the results above listed we suggest that the Casimr vacuum energy can 
be the "glue" that holds together the DNA strands as aggregates. 
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